Effect of umklapp scattering on the magnetic-field-induced spin-density waves in 

quasi-one-dimensional organic conductors 
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We study the effect of umklapp scattering on the magnetic-field-induced spin-density-wave 
(FISDW) phases which are experimentally observed in the quasi-one-dimensional organic conduc- 
tors of the Bechgaard salts family. Within the framework of the quantized nesting model, we show 
that the transition temperature is determined by a modified Stoner criterion which includes the 
effect of umklapp scattering. We determine the SDW polarization (linear or circular) by analyzing 
the Ginzburg-Landau expansion of the free energy. We also study how umklapp processes modify 
the quantum Hall effect (QHE) and the spectrum of the FISDW phases. We find that umklapp 
scattering stabilizes phases which exhibit a sign reversal of the QHE, as experimentally observed 
in the Bechgaard salts. These "negative" phases are characterized by the simultaneous existence of 
two SDWs with comparable amplitudes. As the umklapp scattering strength increases, they may 
become helicoidal (circularly polarized SDWs). The QHE vanishes in the helicoidal phases, but a 
magnetoelectric effect appears. These two characteristic properties may be utilized to detect the 
magnetic-field-induced helicoidal SDW phases experimentally. 



I. INTRODUCTION 

The organic conductors of the Bechgaard salts family 
(TMTSF) 2 X (where TMTSF stands for tetramethylte- 
traselenafulvalene) exhibit a rich phase diagram when 
temperature, magnetic field, or pressure are varied. In 
three members of this family (X=C1C>4, PFg, Re04), 
a moderate magnetic field above several Tesla destroys 
the metallic phase and induces a series o£ SDW phases 
separated by first-order phase transitions.El Because of a 
strong quasi-one-dimensional anisotropy (the typical ra- 
tio of the electron transfer integrals in the three crys- 
tal directions is t a : t b : t c = 3000 : 300 : 10 K), the 
Fermi surfaces of these materials are open. According 
to the so-called quantized nesting model (QNM)Jd^Ll the 
formation of the FISDW phases results from an inter- 
play between the nesting properties of the Fermi sur- 
face and the quantization of the electronic orbits in mag- 
netic field. The wave vector of a FISDW adjusts itself to 
the magnetic field so that unpaired electrons completely 
fill an integer number of Landau levels, thus the Hall 
effect is quantized.™! The standard QNM predicts the 
Hall plateaus of the same sign, referred to as positive by 
convention, which agrees with most experiments. How- 
ever, at certain pressures, a negative Hall effect is also 
observed.E£ry 

We have recently shown that umklapp processes may 
naturally explain the sign reversals of the QHE .113 Our ex- 
planation differs from the one suggested by Zanchi and 
MontambauxMI in invoking the pressure dependence of 
umklapp scattering rather than the electron band struc- 
ture. Although both explanations lead to similar phase 
diagrams, we predict the existence of two SDWs (with 
comparable amplitudes) in the negative phases. More- 
over, we have shown that the negative phases are likely to 
become helicoidal (circularly polarized SDWs) under the 
effect of stronger umklapp scattering. Experimentally, 



this can be achieved by decreasing the applied pressure. 
The helicoidal phases are characterized by, a vanishing 
QHE and a kinetic magnetoelectric effect.EJj 

In this paper we study the effect of umklapp processes 
on the FISDW phases within the framework of the QNM. 
We discuss in more detail the results reported in Ref. [ll], 
and address issues not considered in the latter. The effect 
of umklapp scattering on the FISDW phases was studied 
before by Lebed' using rather crude approximations, but 
the helicoidal phases. &ad the sign reversals of the QHE 
were not discussed. 

In Bechgaard salts, complete charge transfer from the 
molecules TMTSF to the anions X leads to a conduction 
band that is quarter- filled in terms of holes. A dimcr- 
ization along the x axis induces a gap in the electronic 
spectrum. This results in a half-filled band for the holes 
so that umklapp processes transfering Akp — lixja are 
possible (kp being the Fermi wave vector of the holes, 
and a the lattice spacing along the chains). Therefore, a 
quasi- ID g-ology description of the FISDW phases should 
include not only forward (g 2 ) and backward (<7i) scatter- 
ing amplitudes, but also umklapp scattering amplitude 
(53). Since the dimerization is weak, we expect the 
umklapp scattering amplitude (73 to be small. Neverthe- 
less, we shall show in this paper that very weak umklapp 
processes can have drastic effects on the low-temperature 
phase diagram. 

In the next section we obtain the FISDW transition 
temperature in the random phase approximation (RPA) 
and discuss the phase diagram. In absence of umklapp 
scattering, the instability of the metallic phase corre- 
sponds to the formation of a SDW with a quantized lon- 
gitudinal wave vector Qx = 2kp + NG (N integer) 
where G = eHb/h is a magnetic wave vector (with H 
the magnetic field, b the interchain spacing, and — e the 
electron charge). As the field varies, the value of the inte- 
ger N changes (its sign remaining the same), which leads 
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to a cascade of FISDW phases separated by first order 
transitions. The integer N also determines the quantum 
Hall conductivity: a xy = —2Ne 2 /h per one layer of the 
TMTSF molecules.El In presence of umklapp scattering, 

two SDWs with quantized longitudinal wave vectors Q X N ^ 
and Q X ~ N ^ form simultaneously. We label each FISDW 
phase by the integer N such that the SDW with wave 
vector Q X N ^ has the largest amplitude. We find that the 
transition temperature Tc N ^ is determined by a modified 
Stoner criterion that includes the effect of umklapp scat- 
tering. We calculate Tc N ^ numerically as a function of 
g 3 and the detailed geometry of the Fermi surface. In 
order to keep our discussion of the phase diagram com- 
pact, we use some results that are proved in subsequent 
sections. We find that weak umklapp scattering ((73/52 
a few percents) can lead to a FISDW cascade with both 
positive and negative values of N. Since the quantum 
Hall conductivity is still determined by the integer N, 
a xy = —2Ne 2 /h (section VI), this leads to sign rever- 
sals of the QHE as the magnetic field varies. The nega- 
tive phases (i.e., with a sign reversed QHE) correspond 
to even integers N in agreement with experiments. For 
stronger values of 53, the negative phases may become 
helicoidal (i.e., with a circular polarization of the SDWs) 
(section IV). The helicoidal phases are characterized by a 
vanishing QH E an d a kinetic magnetoelectric effect (sec- 
tions VI and VII). There may also be some reentrances 
of the phase N = within the cascade. In section [v|, 
we study how umklapp scattering affects the excitation 
spectrum in the FISDW phases. In section VIII, we con- 



sider the possible coexistence of two successive FISDW 
phases, an issue previously considered by Lebed' in the 
case of sinusoidal SDWs .113 

The experimental cons eque nces of our work are mainly 
discussed in sections O, |VIl| and IX. To a large extent, 



these sections can be read independly of the rest of the 
paper. 



II. INSTABILITY OF THE METALLIC PHASE 
AND PHASE DIAGRAM 

In this section, we first consider the system with no 
electron-electron interaction. We obtain the one-particle 
eigenstates in the presence of a uniform magnetic field H 
along the least conducting axis z, and calculate the bare 
susceptibility. Then we take into account the interactions 
(including umklapp processes) and study the formation 
of SDWs in the IIP A. We also discuss the QHE and the 
polarization of the SDWs using results to be proved in 
the next sections. 

In the vicinity of the Fermi energy, the electron disper- 
sion law in the Bechgaard salts is approximated as (we 
take h = ks = 1 throughout the paper, and the Fermi 
energy fi is chosen as the origin of the energies) 



where k x and k y are the electron momenta along and 
acro ss the one-dimensional chains of TMTSF. In Eq. 
( |2.lD , the longitudinal electron dispersion is linearized 
in k x in the vicinity of the two one-dimensional Fermi 
points ±fcp, and vp — 2at a $m(kpa) is the correspond- 
ing Fermi velocity. The function t±(u), which describes 
the propagation in the transverse direction, is periodic: 
t±(u) — t±(u + 2n). It can be expanded in Fourier series 



t±(u) — —2tbCOs(u) — 2t 2 fcCOs(2u) 
— 2tsb cos(3w) — 2t4b cos(4u) 



(2.2) 



If we retain only the first harmonic (tb), we obtain a 
Fermi surface with perfect nesting at {2kp, tt/6). The 
other harmonics i 2 6, tab " • *C tb generate deviations from 
the perfect nesting. They have been introduced in order 
to keep a realistic description of the Fermi surface despite 
the linearization around i/c^.El In the following, we shall 
retain only tb, t2b, and t^b (as we shall show, t^b does not 
play an important role and can be discarded) . We do not 
consider the electron dispersion along the z axis, because 
it is not important in the following (its main effect is to 
introduce a 3D threshold field below which the FISDW 
cascade is suppressed!!!) . 

The effect of the magnetic field along the z axis is 
taken into account via the Peierls substitution k — > 
— zV — eA. (The charge e is positive since the actual 
carriers are holes.) Following Ref. [2l], we use the gauge 
A = (— Hy, 0, 0). Considering also electron-electron in- 
teractions, we obtain the Hamiltonian Ji = TCq + Ti, ln t 
with 



T-to = ^2 J d2r ^aa{ v )[ v F{-iad x - hp) 



amoj c 



-t±(-ibd y ) + afi B H]ip a(T (r), 



.92 



£ J d 2 ^L(r)V-LW^(r)^(r) 



j2 re -ia>4k F x 



0£,<T 

',t f„\„/,t 



x V'L ( r ) ( r ) (r) (r) • 



(2.3) 



E(k x ,k y ) = v F (\k x \ - k F ) +t±(k y b), 



(2.1) 



Here i/w(r) are fermionic operators for right (a — +) 
and left (a = — ) moving particles, a = +(— ) for f 
(|) spin. We use the notation r = (x,mb) (m integer), 
J cPr = b J dx and a = —a, a = —a. 

Apart from the Zeeman term afisH (fiB is the Bohr 
magneton, and we take the electron gyromagnetic factor 
g equal to 2), the magnetic field introduces the additional 
term arhuj c where rh is the (discrete) position operator 
in the y direction. oj c = Gvf = eHbvp. 

The interacting part of the Hamiltonian contains two 
terms corresponding to forward (32) and umklapp ((73) 
scattering. For repulsive interactions (72,53 > 0. We do 
not consider backward scattering (51), since it does not 
play any role in the mean field theory of the FISDW 
phases. 
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A. Bare susceptibility 

The one-particle eigenstates of H.q were obtained in 
Ref. [H]in the particular case where t±(u) = — 2tbCos(u). 
The extension to a general t± (u) is straightforward, and 
we only quote the final result. The eigenstates and the 
spectrum can be written as 



«„,(>•) = 



./ 1 — m 



f 



du 
27 



I— mi 



-i(l — m)u-\-i-^-T± (u) 



e fc x ,i,<r = v F (ak x - k F ) + aluo c + au B H, 



where we have introduced 



T±(u) = / du't x (u'). 
Jo 



(2.4) 



(2.5) 



L x is the length of the system in the x direction. k x 
is the eigenvalue of the operator —id x (which commutes 
with TLo)- The wave function ; is localized around 
the Zth chain with a spatial extension in the y direc- 
tion of the order of btb/uj c (assuming tb 3> t2b,tzb- ■ ■)■ 
btb/uJc also corresponds to the amplitude of the semi- 
classical orbits in the transverse direction. The localiza- 
tion of the wave functions can be interpreted as Bloch 
oscillations of the electrons in the magnetic field, and 
within the same picture the quantize d s pectrum e£ ; a 
can be seen as a Wannier-Stark ladder £ll The latter pro- 
vides a very natural picture of the quantized nesting 
mechanism, which is-.-a.t .the origin of the FISDWs in 
quasi-lD conductors Indeed, for two rods l\ and 
I2 of the Wannier-Stark ladder, the "nesting condition" 
e fc,,Zi,T = - e t+Q*,hd is fulfilled if Q x = 2k F + (l 1 -l 2 )G. 
Therefore we expect the formation of a SDW at a quan- 
tized wave vector Q x = 2kp+NG (N integer) in presence 
of repulsive electron-electron interactions. 

Now we introduce the bare transverse spin susceptibil- 
ity in the Matsubara formalism 

Xi 0) (r~r',r) = 

(T T 4 T (r,T)Vai(r,r)^ i (r' J 0)^(^,0)) (2.6) 

which is to be calculated with Tlo only. Here r is an 
imaginary time. In Fourier space, we obtain 



xi 0) (q,^) 



1 r r 1 ^ 

1 d 2 rd 2 r' / dre -*i-(«-0+**.r 



L X Ly 



(r-r') 



x X i°)(r-r',r) 
= --^— V f d 2 rd 2 r'e-^ 

L X Ly ^ J 

x G & i ( r ', r , ^G^i (r, r', cj + u„), (2.7) 

where L x L y is the area of the system and q x ~ a2k F . 
G aa is the one-particle Green's function, uj = 7rT(2n + l) 



(n integer) is a fermionic Matsubara frequency, and 
uj v = 2-kTv [y integer) a bosonic Matsubara frequency. 
Using 

C,(( r )C/( r ') 



G aa (r,r',u) = 



~As~.J,i 



we obtain the well known resultS 



(2. 



(2.9) 



where x a L '(fei^) is the susceptibility of a one- 
dimensional system without interaction. In the static 
limit (uj u = 0), 



X a D (lx) 



N(0) 



In 



(2jEo 
1 ttT 



-Re* 



v F 



(q x ~ a2k F ) 



(2.10) 



where N(Q) — 1/irvpb is the density of states per spin, 
* the digamma function, and 7 ~ 1.781 the exponen- 
tial of the Euler constant. Eo is an ultraviolet cut- 
off of the order of the bandwidth. Since x^i^kp) = 
(N(0)/2)ln(2~fE a /TrT) 1 the bare susceptibility x (0) has 
logarithmic divergences at quantized values Q x = 
2kp + NG (N integer) of the longitudinal wave vector. 
The coefficients I n {lv) defined by 



— imqybr— r+ 

J 1 — mJ V —m 



2tt 



du 
— < 
2tt 



l(l-l')u+ T ^[T ± (u+q v b/2)+T ± (u-q y b/2)} 

(2.11) 

are well known in the QNmE They crucially depend on 
the detailed structure of the quasi-lD Fermi surface and 
therefore determine the stability of the metallic phase 
with respect to the formation of a SDW phase. 



B. RPA susceptibility and phase diagram 

We now consider the total Hamiltonian as given by 
fl2.3p , and introduce the spin susceptibility 

X aa >(r,r',T) = 

(T T 4 T (r, r)^i(r, r)^, x (r', 0)^ T (r', 0)). (2.12) 

In the RPA, Xaa' is determined by the integral equation 
(see Fig. |) 

w(ry) = Mf (r-r') 

+.92 J d 2 r 1 x { ° ) {r~r 1 ) Xaa '{rur') 

+53 / d 2 rix£Hr-riy a4kFXl Xaa'(ri,r') 

(2.13) 
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where we now consider only the static limit (lo u = 0) . In 
Fourier space, this leads to 

+g2X { a ) (Qx)Xaa'{qx,q' x ) 

+93X i " ) (qcc)x& a '(q x ~ aAk F ,q' x ). (2.14) 

We have not written explicitly the dependence on q y since 
the latter is a conserved quantity. Because loga- 
rithmically diverges at Qx = 2kp + NG, we consider 
spin fluctuations only at these wave vectors. Because of 
umklapp processes, fluctuations at Qx are coupled with 



fluctuations at Qx 



(N) n (ATh 



4fc F 



x ++ (Q x N \Q 



xf(Qi N) ) 



x-+(-Q 



92xf{Qi N) )x ++ {Qi N \Qi N) ) 
g,xf(Q { x N) )x- + (-Q x N \Q x N) l 
" r>'"^= 92X ^\-Q^) X . + (-Qf\QW) 

W(-Qf)) x++ (QW,QW). 

(2.15) 



+93X- 



Using X&\~qx,q y ) = Xa'(qx,q y ), we obtain 



x++{Qn,Q 



N 



xT(Qn)[i 



-92X { +\Qn)} 



x-+(-Q { /\Q v ;Qn) 



[l-92xf (Qn)][1-92X { + ) (Qn)} 
„ .io),~ (o), 



■ 5ix+ 0) (Qw)xi 0) (Qjv)' 



93X ( + ) (Qn)x ( + ) (Qn) 



[l-52xf(Qiv)][l-52xi 0) (Qiv)] 



9zX+ (Hn)X+ 



Qn) 



(2.16) 



We have written explicitly the dependence on the 
transverse wave vector by introducing 



Qa^ = (Q ( x \Qy) and Qjy - (QT^-Qy). (2.17) 

Note that in our notation the wave vector Qjy = Q n 
has both the signs of N and Q y reversed compared to 
Qat- This happens because Umklapp scattering couples 



Qx ,Qy), but the latter is equivalent 

to (Q x N \~Q y ). 

In presence of umklapp processes, the transition tern- 



(Qx N) ,Q v ) to (- 



perature T< 
terion 



(AT) 



is determined by the modified Stoner cri- 



1-92X ( +\Qn)][1-92X [ + > (Qn)} 



(0), 



g!x i + ) (QN)x y :'(QN) = o, 



(0), 



(2.18) 



which is the condition for vanishing of denom inators and 
divergence of susceptibilities Xaa' in ( [2.16 ). Two fluc- 
tuation modes diverge simultaneously, which leads to 
the formation of two SDWs with wave vectors Qat = 

(Qx N \Q y ) and = (Q { x N) ,-Q y ) (see (p^) ). We 
label each FISDW phase by the integer N such that 
the SDW with wave vector Qat has the largest ampli- 
tude. This is equivalent to defining N by the condition 



X 



(0) 



(Qn)>X+\Qn) (section 0. 



We will show in section |v| that the SDWs can be ei- 
ther sinusoidal or helicoidal, depending on the value of 
the angle tf N (e] - 7r/4, 7r/4]) defined by tan(2#/v) = 
2B/(A N - A N ) and 



.4 



1 



±N 



T 2 

X ±N 



92 



B = — 



(Qy) \9l-9l 
93 



X ( + ) (Q±n)\, 



In{Q v )In(Qv){92- S3) 



(2.19) 



r 



For sin 2 (2??Ar) < 2/3, the SDWs are sinusoidal and the 
Hall effect is quantized: a xy = —2Ne 2 /h. N then cor- 
responds to the quantum number which is directly mea- 
sured in transport experiments. For sui 2 (2i9at) > 2/3, 
the phase is helicoidal. The Hall effect vanishes, but a 
kinetic magnctoelectric effect appears (see sections VI 
and VII). The phase N = is sinusoidal if Q y = ir/b 
and helicoidal if Q y 7^ n/b (section IV). -On also deter- 
mines the ratio I7I (0 < I7I < 1) of the amplitudes of the 
two SDWs: 



-fjy(Q a )tan(#Af) - rIfj(Q y ) 



(2.20) 



where r = 53/(72- The figures of this section also show 
I7I = I tan($jv)| which is analogous to I7I but for the ef- 
fective mean-field potential acting on the electrons. (7 
and 7 are precisely defined in section IV.) I7I and I7I 
increase with |$aH- At the transition between the sinu- 
soidal and helicoidal phases I7I ~ \j\ ~ 0.518. 

In the absence of umklapp processes, Eq. ( [2.18 ) yields 
the usual Stoner criterion 1 — (Qn) = for the for- 

mation of a SDW at wave vector Qat = (Q X N \ Qy)- The 
quantized longitudinal wave vector Qx = 2kp + NG 
and the transverse wave vector Q y are chosen to max- 
imize the transition temperature Tc at a given mag- 
netic field. Except when N — 0, Q y is incommensurate: 
Q y ^ ixjb. The SDW is sinusoidal (7% = for g 3 = 0), 
and the quantum Hall conductivity in the FISDW is de- 
termined by the integer N: a xy = —2Ne 2 /h. As the 
magnetic field increases, the value of N changes, which 
leads to a cascade of FISDW phases separated by first 
order transitions. In the simplest version of the QNM, 
where t^b = t±b — ■ • ■ = 0, the phases of the cas- 
cade are labeled by N — •••5,4,3,2,1,0 as the field 
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increases. The integer TV is positive provided t 2 b > 0, 
i.e., sgn(TV)=sgn(i 2 b). (For t 2 b < 0, one would obtain a 
similar sequence but with negative values of TV.) 

We study the phase diagram in presence of umklapp 
scattering numerically (see Figs. §-?])■ The calculations 
are done for tb = 300 K, t 2 b — 20 K, t^b = K, and 
E = 2000 K. g 2 + g 3 = 2/ \n(2-fE /nT^) is held fixed 
where T c °° = 12 K is the transition temperature for an 
infinite magnetic field, and gt = N(0)gi (i = 2, 3) are di- 
mensionless coupling constants. The figures are obtained 
for different values of t 4 b and the ratio r = gz/g 2 . The 
parameters are such that we are in the weak coupling 
limit: g% < g 2 < 0.4. In Bechgaard salts, t 4 b is expected 
to be-i a very small energy scale, of order 1 K or even 
lessB 

For r — and t 4 b = 0.75 K, we obtain the sequence 
N = •••,4, 3, 2, 1,0 in agreement with what has been 
found for t 4 b — (i.e., a small value of t 4 b does not change 
the phase diagram when r = 0).EI The transverse wave 
vector Q y varies approximately linearly with the field 
within each phase, and is incommensurate (Q y 7^ tt/6) 
except in the phase N = (Fig. ||). 

A very small 53 does not change the phase diagram 

qualitatively compared to the case 33 = 0. Now the 

main SDW at the wave vector Qjv coexists with a weak 

SDW at the wave vector Qjy. In general, the value of Q y 

that maximizes X (Qiv) does not maximize X (Qjv)i 

so X (O) (Qa0 < X (0) (Qjv)- As a result ' the SDW am- 
plitude at the wave vector Q^y is very small, and the 

polarizations of the SDWs are linear. The values of N 
follow the usual positive sequence N = • • • , 5, 4, 3, 2, 1, 
as the magnetic field increases. 

A larger value of 33 increases the coupling between the 
two SDWs. This leads to a strong decrease of the criti- 
cal temperature or even the disappearance of the SDWs. 
However, for even A, there exists a critical value of 53 
above which the system prefers to choose the transversely 
commensurate wave vector Q y = n/b for both SDWs. 
The reason is that, for even N (as opposed to odd N), 
Q y = n/b corresponds to a local maximum of the suscep- 



tibilities and X {Q) (Qi N \n/b) ~ x m (Qx"',Tr/b). The two 
susceptibilities are strictly equal at to, = 0, but when 
Ub > 0, X {0 HQ ( J barN \*/b) > X {0 HQi N) ^/b) (this re- 
suit also holds fpr t 3 b =/= since x^ {Qx , t/6) is inde- 
pendent of t3b)£3 This yields a negative Hall plateau, 
provided the SDWs are sinusoidal. Thus, for r = 0.025 
(g 2 ~ 0.37 and gz — 0.01) and t 4 b = 0.75 K, we find 
the sequence N = ■ ■ ■ 5, 4, 3, -2, 2, 1, (Fig. g). A neg- 
ative commensurate phase with N = — 2 and Q y = n/b 
appears in the cascade. All the phases are sinusoidal, 
so the Hall effect is quantized (a xy = —2Ne 2 /h). For 
r = 0.025 and t&b = —0.75 K, we obtain only positive 
integers: N = •••5,4,3,2,1,0 (Fig. |). However, the 
phase N ~ 2 has split into two subphases: for large 
enough field, the transverse wave vector Q y differs from 
7r /b and varies linearly with the field; for a weaker field, 
the phase is commensurate (Q y = n/b). All the phases 



are sinusoidal (sin 2 (2i?Ar) < 2/3) so that the Hall effect 
is quantized. 

Fig. |H shows that both SDWs have comparable am- 
plitudes in the negative phases (and more generally 
in the phases where Q y = n/b (Fig. ||)): |7|,|7| — 

0.3. This results from the property x^°\Qx^/b) — 



X 



^{Qx^^n/b). On the contrary, the amplitude of the 



SDW with wave vector Q^y remains very small in the 
positive phases. Note that it is 7 which is actually van- 
ishingly small, \ j\ being of the order of r (see section IV A 
for a further discussion). 

The strength of umklapp scattering is very sensitive to 
pressure. Indeed, hydrostatic pressure reduces the dimcr- 
ization gap and diminishes 33. Therefore, we conclude 
that sign reversals of the QHE can be induced by vary- 
ing pressure. In our simplified model, this effect requires 
tib > 0. Our results provide a new explanation of the sign 
reversals of the QHE which have,long been observed in 
quasi- ID organic conductors Jl3~tJ In particular, Balicas 
et al. have recently shown unambiguously the existence 
of the phase N = -2 in (TMTSF) 2 PF 6 at a pressure of 
8.3 kbar by observing a sign reversal of the QHE with, 
a well defined Hall plateau corresponding to N — — 2. til 
(These results required a conditioning procedure in which 
current pulses depin the FISDW from lattice defects and 
tend to reduce hysteresis.) The observed FISDW cascade 
corresponds to N = ■ ■ ■ 4, 3, —2, 2, 1, 0. When the pres- 
sure is increased to 9 kbar (which decreases the umklapp 
scattering strength), the phase N — —2 disappears and 
the usual sequence N — ...4, 3, 2, 1, is obtained. 

If the value of t±b is reduced, the phase N — —2 be- 
comes helicoidal. This is shown in Fig. || obtained for 
r = 0.03 and t 4b = 0.3 K. (When t 4b = 0.3 K, the 
phase N = — 2 appears for stronger umklapp scatter- 
ing. This is the reason why we show the phase dia- 
gram for r = 0.03 and not r — 0.025.) In the he- 
licoidal phase, I7I ~ I7I > 0.5. For t 4 b = 0, there 
is a degeneracy between N and — N at Q y — n/b: 

<°>(QW */6) = xf (QfW&), iN^/b) = I 



N^/b) 



X 

and Am = Aff. This yields i?jy = 7r/4 and the waves 
are helicoidal. A finite t 4 b lifts this degeneracy, so that 
A n ^= A^ and in turn sin 2 (2i9jv) < 1. Thus, the stability 
of the helicoidal phases is strongly related to the degen- 
eracy between TV and — N occuring for Q y = n/b. In our 
model, this degeneracy is entirely controlled by t 4 b- It 
is not affected by t%b, which is the reason why we have 
chosen t^b = in the numerical calculations. 

If, on the other hand, the value of t 4 b is increased, 
the ratio of the amplitudes of the two SDWs decreases. 
For instance, for r — 0.025 and t 4 b = 1.5 K, we find 
|7|, |7| ~ 0.15. Thus, for large t 4b (i.e., t 4b > 1.5 K), the 
amplitude of the SDW with wave vector Q^y becomes 
very small. There is then no real difference between a 
positive and a negative phase (except for the sign of the 
QHE) insofar as both contain a main SDW, which coex- 
ists with another SDW with a very small amplitude. 
If r is increased to 0.06 (with t 4 b — 0.75 K), a sec- 
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ond negative phase (N = —4) appears, and the cas- 
cade becomes N = • • • 8, 7, -4, 6, 5, 4, -2, 2, 1, (Fig. §). 
Note that N = —2 and iV = — 4 are the two, nega- 
tive phases most easily observed in experiments. t3[3 The 
phase N = — 2 has split into two subphases: one is he- 
licoidal (sin 2 (2tf 2 ) > 2/3), one is sinusoidal (sin 2 (2^2) < 
2/3). Thus, increasing the strength of umklapp processes 
makes the negative phase N = — 2 helicoidal. 

In order to observe the helicoidal phase experimen- 
tally, it would be desirable to stabilize the negative phase 
N = — 2 at the lowest possible pressure (which corre- 
sponds to the strongest g 3 ). In (TMTSF) 2 PF 6 , the pres- 
sure has to be larger than 6 kbar, since below this pres- 
sure the FISDW cascade disappears^! In the experiment 
reported in Ref. [l4|, where the phase N — — 2 has been 
observed at 8.3 kbar, the pressure could be reduced only 
by about 2 kbar. Nevertheless, because 53 is very sensi- 
tive to pressure, such a pressure reduction could induce 
a significant increase of the umklapp scattering strength. 
(TMTSF) 2 Re0 4 , where sign reversals of the QHE have 
been observed under pressure,E3 could also be a good 
candidate for the observation of helicoidal phases. In 
(TMTSF) 2 C10 4 , sign reversals-of the QHE have been 
observed at ambient pressure, Iffil so that it is not possi- 
ble to increase 53 by decreasing pressure. 

When r is further increased, only phases with nega- 
tive even N survive. This leads to the sequence N — 
■ ■ ■ — 8,-6,-4,-2,0 for t^, — 0.75 K apart from some 
reentrances of the phase N = within the cascade to 
be discussed below. (For t^b = —0.75 K, we obtain the 
sequence N — •••8,6,4,2,0.) Moreover, all the phases 
N ^ are commensurate (Q y = n/b). This is shown in 
Fig. @ obtained for r = 0.2, i.e., g 2 ^ 0.32 and 53 ^ 0.06. 



The phase N = is somehow special since X+ (Qn) 



(N 7^ 0) which decreases 
n/b and X f (Q&*° , tt/6) = 



X+ (Qjv) m that case. As a result, the transition tem- 
perature ri°^ is determined by 1 — (p 2 + 93)x+ ^ (Qo) = 
and does not depend on the ratio r = g-i/gi when 
92 + 33 is held fixed. This should be contrasted with the 
transition temperature T c 
with r (except when Q y - 

X+\Qx N \ ft / b)) . This explains why, when r increases, 
some reentrances of the phase N = are observed within 
the cascade. Notice that the latter phases are not com- 
mensurate (Q y ^ n/b) contrary to the last phase (N = 0, 
Q y = n/b) of the cascade. The reentrant phases N = 
are always helicoidal (independently of the structure of 
the Fermi surface), but the last phase N = 0, Q y = n/b 
of the cascade is sinusoidal (section IV). 



C. Effect of ID fluctuations 

Our numerical results show that a very small value of 
the umklapp scattering amplitude is sufficient to explain 
the phase diagram of Bcchgaard salts. In fact, the order 
of magnitude of gs which is required to stabilize negative 



phases strongly depends on the choice of the ultraviolet 
cutoff E . 

Within a mean-field picture, Eq is a large energy, of the 
order of the electron bandwidth. It has been argued that 
because of ID fluctuations, the appropriate cutoff to be 
used in the QNM is not the bandwidth but the dimen- 
sional crossover temperature T x i <C EqL3 (Above T x i, 
the behavior of the system is essentially ID, so that the 
interference between particle-particle and particle-hole 
channels invalidate the mean-field (or ladder) approxima- 
tion.) A stronger value of g$ is then required to stabilize 
negative phases. This is in agreement with the sugges- 
tion of Behnia et al. that the effective low-temperature 
value of (73 in a magnetic, field is significantly enhanced by 
ID high energy scalesBj This point of view is supported 
by NMR measurements and a large magneto-resistance, 
which shows an activated behavior becoming more and 
more pronounced as the field increases.t3 



III. ORDER PARAMETERS: HELICOIDAL VS 
SINUSOIDAL WAVES 

The divergence of the susceptibilities \aa' (section ||) 
indicates that the FISDW phases are characterized by 
the order parameters 

A Qa (r) = (^LWfe(r)) = £ A^e'^-. (3.1) 

P =± 

The two wave vectors Q p n, P = ilj are given by Eq. 



@. Q 



,(pJV) 



2k F + pNG and Q y N) 



-Q 



(AO 



The 



= complex numbers Aq„ are the order parameters of 



the FISDW phase. A Q(T (r) = A^(r) implies A 



(pN) 



A y - F - 



(pN)* 



Among the eight order parameters A„a , only 
four are therefore independent and sufficient to character- 
ize the SDW phase. Note that for N = 0, one should dis- 
tinguish in general between the phases N and N. In this 
case, there are two SDWs with wave vectors (2fcp, Q y ^) 



and (2kp, Q^) with <^ :/ 



(S) _ 



= -Q y Q) . For Q 



(( V 0.7T/6, 



the two SDWs are different so that both Aaa and Aaa 
are needed. When Q y ^ = 0, n/b, the two SDWs are iden- 
tical, and only one order parameter (for instance Aaa ) is 
sufficient. 

Now we discuss how the polarization of the wave affects 
the order parameters A^a^ . For simplicity, we consider 
only one wave vector Q and denote the four different 
order parameters by A aa (among which only two are in- 
dependent since A aS = A* CT ). For a SDW polarized per- 
pendicularly to the magnetic field, the expectation value 
of the spin density operator S(r) can be written as 

a,tr,cr' 

= m x cos(Q - r + 6i), 
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(Sy(r)) = £ (V'L(')^Vaa'(r)) 

= m y cos(Q • r + #2), 



(3.2) 



where we have introduced 

A aa {r) = g 2 A aa (v) + g.e-^^A^ir). (4.2) 



where r (x) and are Pauli matrices. 9\ = 9 2 corre- A Qcr (r) is given by ( |3~T| ). For g 3 7^ 0, the mean-field 
sponds to a sinusoidal wave, while 61 — 9 2 ± 7r/2 and potential A Q(T (r) acting on the electrons is a linear com- 
m x = m y corresponds to a helicoidal wave. (S(r)) can bination of the order parameters A aa (r). 



be expressed in terms of the order parameters ( p.l| ): 
(S+(r)) = (S x (r))+i(S y (r)) 

= £<Vi T (r)^i(r))=£A aT (r). (3.3) 

a a 

Comparing this expression with Eqs. ( |3.2] ), we find 
1 



Calculating the free energy (per surface unit) of the 
phase N to the fourth order in the order parameters, we 



obtain Fn 

p(2) 



F 



(2) 



N 



F^ ] with 



r N — 



E 



d 2 r - 



L X Ly 



A QT (r)A ai (r) 



T 



A +T = - (m.e-^+Jmje-' 92 ) 



L X Ly 



J d 2 r 1 d 2 r 2 K a] {Y 1 )K &l (v 2 ) 



A_| = - (m^e^ 1 + im 9 e lS2 j 



(3.4) 



For a sinusoidal wave |A+j| = |A_j|. For a helicoidal 
wave A_|_| 7^ and A_j = (or the symmetric solu- 
tion A_j 7^ and A+f = 0). The inverse is also true: 
|A+f| = |A_f| implies that the wave is sinusoidal, while 
A+t 7^ and A_| = implies that the wave is helicoidal. 



IV. POLARIZATION OF THE PHASE TV 

In this section we derive the Ginzburg-Landau expan- 
sion of the free energy as a function of the order parame- 
ters At N) - The minimum of the free energy determines 
the polarization (linear or circular) of the SDWs. 

The mean-field (or Hartree-Fock) Hamiltonian is given 

by 

Hmf =n -J2 [ d 2 rA Qa (r)4 ff (r)?/w(r) 

+ £ J d 2 rA QT (r)A al (r), (4.1) 



P (4) _ 
N - 



x G qT (r 1 , r 2 , w)Gai ( r 2, ri , a;), 

T 
2L 



^7~E / rf 2 nd 2 r 2 d 2 r 3 d 2 r4A QT (ri)A gi (r2) 



xA a] -(r 3 )A al (r i )G a] (r 1 ,r 2 ,u})G &l (r 2 ,r 3 ,Lu) 
xG QT (r3,r4,w)G a x(r4,ri,w). (4.3) 

G a a is the s ingl e-particle Green's function in the metallic 
phase (see (^8)). Introducing the electron-hole pairing 
amplitude (we use the notation Wj = (kj X , lj) for the in- 
dices of the functions 4>%_ v (Eq. (2.4))) 



A aa { Wl ,w 2 ) = J d 2 rC*(r)C 2 (r)A QCT (r) 
= A%^(w 2 ,wi), 
we write the free energy as 
d 2 r ~ 



(4.4) 



F 



(2) 



N 



E 



L X Ly 



A aT (r)A ai (r) 



T 



J 



L X Ly 



E 



a,u),wi ,W2 



A aT (wx, w 2 )A al (w 2 , wi) 



F 



(4) 



T 



N 



E E 



A a ^(wi,w 3 )A &l (w 3 ,W4)A a ^(w i ,w 2 )A al (w 2 ,wi) 
( iuJ - e Sa)(^ " " <Lu)(^ - Ct) 

I 



(4.5) 



(4.6) 



sfP = 



b.n7r 



The pairing amplitudes are given by (appendix |A|) 
A Q ,K, W2 ) = £ S kMaQipN) ( 92 A^ + 53 Ar ) ) 

-iaQ^ N) b{h+i 2 )/2 T jr>( N )\ (A 7\ A± N and B are defined by (|2.19|). We have introduced 



p=± 



^(A^+c.c.) 

(4.8) 



the new order parameters A^^' related to Aq/' by 



This leads to 



F 



( P N), 2 



BA 



at 



SF f 



(2) 



' 



At N) = I pN (92At N) + 93 Al N) ) = A^> 

g 2 IpNA a P a^ — g 3 I p NA a ^ 



A^ N) 



(4.9) 
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with I pN = I p N{Q { y N) )- 

A commonly used approximation in the QNM is the 
quantum limit approximation (QLA) valid when lu c 3> 
tMIlI It consists in retaining only the most singular 
(electron-hole) pairing channels that have the logarithmic 
singularity ~ \n(2jEo/irT). This singularity results from 
pairings between electron and hole states of the same en- 
ergy. Therefore, in the QLA, A QCT (u>i, w 2 ) is nonzero only 



if e a - 



This leads to 



scribed with only one gap, although gaps opening above 
and below the Fermi level are still neglected (see section 
^). For this reason, we use the term QLA rather than 
SGA. In the QLA, we would have obtained Eq. (4.S), to- 
gether with ( [2.19 ), but with the exact susceptibility x+ 
replaced by 



xf(Qpiv) 



QLA 



j2 m ln 



(4.11) 



A aa (w 1 ,w 2 ) 



QLA 



UK, 

P =± 



Xe -i a Q^)b(2h- a pN)/2^(pN) _ ( 4 _ 10 ) 



The QLA is usually known as the single gap approxima- 
tion (SGA) because it amounts to considering only the 
gap that opens at the Fermi level. However, when umk- 
lapp processes are present, the spectrum cannot be de- 



J 



Although qualitatively correct, this approximation 
strongly underestimates the susceptibility x^ ' because 
it neglects terms ~ ln(2^fE /\n\uj c ) (n ^ 0) with respect 
to ln(2jEo/irT). However, the QLA becomes very accu- 
rate for higher order contributions to the free energy. For 
Fff, corrections to the QLA are of order T 2 /uj 2 and can 
therefore be neglected when uj c » T. Within the QLA, 
we obtain 



r AT Z 



K 



N 



2 - E 

a Pi-.P2,P3,Pi 



a] 



A (Pa N ) * I ( P3 A) X (p 4 A) * 

at a] a] 



x exp 



x6 (Pi-P2+P3-Pi)N,0 ^gln") 



2p 2 



[iN{b/2)[ Pl Q^ - (2 Pl -p 2 )Q^ + {2 Pl 

Q i y P2N) +Q i y P3N) -Q i y PiN \n2^/b^ 
I 



p 3 )Q { y P3N) 



(P4A)1 



where K = 7C(3)A(0)/(16tt 2 T 2 ) and C(3) ~ 1.20. A 
somewhat lengthly calculation (see appendix III) leads to 



F 



(4) 



K 



A 



= tEK! 



(pJV),4 



SF 



(4) 



2 

+A (0) 



i \ 1 v(iV)|2 



P ■ n 



b.mt 



a 

I 



SF 



(4) 



«1 



+ 2 (A -r 



'at 



) 2 +C.C. 



(4.13) 



(AO 
"1 

at 



COs(z?jv) 

sm(d N ) 



- sin(i?jy) 

COs(^at) 




we obtain 



(4.12) 



(4.14) 



4 2) = E(^7 



A^* 



.4 



A 



E 



(AT) 

«I 

(A) 
at 



{\[ N) \u^ + \ 2 N U N) \ 2 } (4-15) 



Eqs. ( |4.8| ) and (|4.13[ ) show that the case A = and 



)! ' — 0,7r/6 is special, since in that case one cannot 



distinguish between and A^ 



In practice, the case 

= never occurs (see section [0]) so that we shall 
not consider it any more. 



for 



tan(2$Ar) 



2B 



A 



N 



A 



A 



(4.16) 



2$jv can be chosen in ] — n/2, tt/2]. The eigenvalues \\ 



(A) 



and are given by 



A. Phases N / 0, or N = and 



We first consider the cases N ^ 0, and N — with 
ir/b. The latter corresponds to the reentrant 
phases N — within the cascade ( see section |l|) . The 
quadratic part of the free energy (L8) is not diagonal 



in the order parameters since A^ is coupled to A^ 
:ing the new order parai 

fined by 



(A) 



Introducing the new order parameters and de- 



A N cos 2 ($ n ) + A N sm 2 ($ N ) + Bsin(2'd N ) 
A N + A N 1 



rSgn(Ajy - A N ) 



A^f + AB 2 ] 1 / 2 , 



2 

x[{A N 

A N sin 2 (i?Ar) + A N cos 2 (#jv) - Bsm(2tf N ) 
A N + A N 1 

2 2 Sgn< - N ~ ^ 

x[(A Ar -^) 2 + 4B 2 ] 1 / 2 . (4.17) 
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The transition temperature Tc N ^ is determined by 
min(A^ Ar \ A^) = 0. In this section, the sign of TV has 
not been specified yet, since both phases iV and — TV 
correspond to a phase with two SDWs at wave vectors 

(Qi N) ,Qv N) ) and (Q { f\-Q ( y N) ). In the following, we 
fix the sign of TV by the condition x ( +(Qn) > X+\Qn) 
(i.e., in the absence of umklapp processes, the phase TV 
would be more stable than the phase —TV). Then we have 



< A N < Aft and < \{ N) < A^ JVj for T > Tk"> (as- 
suming < p3 < gi)- The transition temperature T ( 



(AO 



(N) 



is determined by X{ 



(N) 



0, i.e., A N A n = B 2 . Using 



Eqs. ( pl9| ). we recover the result ( |2.18| ) obtained from 
the RPA calculation of the susceptibility. 



- X ( 2 N) > for T < T C (A °, we can put v { a N) = 

the expansion of the free energy. This leads to 



K 



COS 4 (^at) +sin 4 (^Ar)] 



x \uW\ 4 + Ksm 2 (2{) N )\u^ ) u {N) \ 2 . (4.18) 

The minimum of the free energy is determined by 
<>F " .W \ X ™ + X[cos 4 (^) + sin 4 (^)]| U W| 2 



For \{ N) < (i.e., T < T t 
corresponding to a SDW phase: 

i) | U W| = (^ 0). Using ,r 



+i^sin 2 (2 l 9 A r)|4 Ar) | 2 l =0 



(4.19) 



I AJ? | = 



| AW | 



there are two solutions 

0, we find 
and also |A^| = 



|AW| = 



A^|. Fr om the results of section tt IU 
we conclude that the two SDWs are linearly polarized. 

(or the symmetric solu- 



ii) «w f o and u w 



tion: u_ 7^ and — 0). Va — then implies 



A^ = A^ - 0, and A^ = 



A ( ^ = 0. 



This corre- 
sponds to two helicoidal SDWs of opposite chiralities. 
For the sinusoidal waves (solution (i)), we find 



,W| 2 = — 



(A') 



X I 



±sin 2 



(20*)) ' 



(4.20) 



A 



X(l + isin 2 (2^Ar))' 
while for the helicoidal waves (solution (ii)), we have 



(4.21) 



,W|2 



(JV) 



if 1 



isin 2 (2^)) 
,(W) 2 



u W = 0, 



AT 



2K (1 



! (2tf7v)) 



(4.22) 



(4.23) 



Comparing Eqs. (4.21) and (4.23), we conclude that the 
helicoidal phase is more stable than the sinusoidal phase 



when sin (2i?jv) > 2/3, i.e when y/2\B\ > \A N - Aft\. 
The free energy Fn (Eq. ( 4.18 )) is shown in Fig. || as 

a function of and \u_ | for A^ < 0. When 

sin 2 (2^Ar) < 2/3, the minimum of Fjy corresponds to 
When sin 2 (2$ at) > 2/3, there are two 



minima located on the lines u+ J = and = 0. For 
sin 2 (2^Ar) — 2/3, the minima are infinitively degenerate. 

In the absence of umklapp processes (g^ = 0), ^at = 
and the SDWs are sinusoidal. For values of g% such that 
Qy — ir/b, the polarization of the waves depends on 
the details of the Fermi surface. In our model, it is de- 
termined by tib (see the numerical calculation of section 

The results obtained in this section are very simple 



when TV = 0. Since X 1 +(Qn) = X+'(Qn) for N = 



(o), 



(the susceptibility does not change when Q y 



An = Aft and $ = tt/4. The phase TV = (with 
Q y ^ 7r/6) is therefore always helicoidal. This result 
does not hold when Q y = ir/b (see section IV B). 

Using — 0, we obtain the ratio of the amplitudes 
of the two SDWs: 



aw 



tan($jv), 



jI N - rl 



N 



In 



rjI N 



(4.24) 
(4.25) 



Since d N e] - 7r/4,7r/4], | 7 |,| 7 | < 1. The SDW with 
wave vector Qat has the largest amplitude. This sim- 
ply results from the fact that we have chosen TV such 
that x (0) (Qat) > x (0) (QaO- The transition from the 
sinusoidal phase to the helicoidal phase occurs when 
sin 2 (27?A,) = 2/3, i.e., when \^\ = 0.518. 

Note that in the positive phases, it is 7 and not 7 which 



is vanishingly small. When X+\Qn) ^ X+'(Qn), the 
minimum of the free energy corresponds to A 



(AT) 







0), since it is the potential A a | (and 



(and not A 

not A Q |) that couples to the electrons (see Eq. ([4.l|)). 
7 ~ implies | 7 | ~ r. 



B. Phase TV 



0, Q< 0) 



: n/b 



We now consider the case TV = with Q y 0) = n/b, 
which corresponds to the phase TV = terminating the 
cascade. There is only one SDW in that case, so we 
can choose AqJ = A«j (alternatively, we could choose 
Ai°2 = 0). The free energy is then given by 



Fn =4 



EhlAi?! 2 



BA^Af}* 



2K\ A^ | 4 



(4.26) 







The quadratic part of the free energy Ffi is diagonalized 
by introducing the order parameters u and v defined by 



(4.27) 




This leads to 



with 



p(2) 



x (0) _ 1 



4Al°V| 2 



(4.28) 



1 



A 



U \92 + .93 
,0) = 1 ( 1 



-xf(Qo) 
-xf(Qo) 



(4.29) 



Tr°^ is determined by Aj° 



0, i.e., 1 - ( ff 2 + 
ff3)x+ 0) (Qo) = 0. This agrees with the RPA result ( pTS| ) 
when x+^(Q/ 



xf (Qiv) = x+ j; (Qo) 



Since A 



(0) 



(0), 

for T < T c (0) . This implies 
— iimo. uuc jjxiase N — 0, Q y — ir/b is always 

linearly polarized. For N = (and Qy '' = 7r/6), g?, cou- 



> A^, » 



A (0 -> = AL ]. Thus, the phase N 



pies A^ } and A^. 



As a result, we cannot have A^°j ^ 



and A^°] = 0, so that the polarization cannot be circu- 



lar. The case g% = is special since then = \£ J and 
both u and v become nonzero below T c . Nevertheless, 



,(0) 



Eq. (4.26) shows that the sinusoidal wave is more stable 
than the hclicoidal wave in this case too. 
Setting v = 0, we obtain the free energy 



^o 



4Al 0) M 2 



AK\u\ 4 . 



Minimizing Fq with respect to u, we obtain 



A^ 
2K 



Fn = - 



(aT)! 

K 



(4.30) 



(4.31) 



V. SPECTRUM 

In the absence of umklapp processes, the spectrum con- 
sists of Landau subbands separated by gaps opening at 
the Fermi level, and,.at-jria; c /2 (n integer) above and be- 
low the Fermi leveLBETLl We show in this section that for 
sinusoidal SDWs the gaps opening at the Fermi level de- 
pend on the transverse momentum k y (in this section, we 
use a gauge where k y is a good quantum number). On 
the other hand, the spectrum of the helicoidal phase re- 
mains independent of k y . As shown below, this provides 
a natural explanation of the instability of the sinusoidal 
phase with respect to the helicoidal phase. 

In this section we use the gauge A = (0, Hx, 0) where 
the calculations are simpler. In the next section this will 
also allow us to calculate the QHE using previous results 
derived in this gauge. Using the Peierls substitution, we 
obtain 

Ho = ^2 J dx^ acr (x,ky)[v F (-iad x - k F ) 

a,cr,ky 

+t^(k y b - Gx) + a^ B H]ip aa (x, k v ), (5.1) 

where we have taken the Fourier transform with respect 
to the y direction. In the absence of electron-electron 
interaction, the eigenstates and the spectrum are given 
by 



1 

7F 



ik-r+i-^-T ± (k y b-Gx) 



e k cr = v F (ak x — k F ) + <j\i B H. 



(5.2) 



The dispersion law is now one-dimensional and the states 
are extended in both the x and y directions. However, 
since the energy does not depend on k y , we can take the 
Fourier transform with respect to k y and obtain local- 
ized wave functions similar to those obtained in section 
0. The Wannier-Stark ladder can then be recovered by 
a simple gaug e transformation. 

Using (4.1), the interacting part of the mean-field 
Hamiltonian can be written as 



J 



U 



MF 



„ 7, J 



V (JV) 

j—'ipt l (x,k y )<ip a i(x,k y + aQ y N) ) 



{N)\ 



AN) 



- A K 



■f^-i>L(x,k y )'il>cii(x,k y +aQ y N) ) 

N 



h.C. 



(5.3) 



up to a constant term. Introducing new fields defined by 



ih (v h \ ~ P iakFX ~ iaa!L ^ x+i ^ T ± i - k y b ~ Gx \l ( <t» b \ 



(5.4) 



we rewrite the Hamiltonian as 



^mf= ^2 / dxip f Q(J (x,ky)(-iav F d x )?p a(7 (x,ky) - ^ j dx 



j~iuNGx-ia7i{k y b~Gx+aQ[ N) b/2) 



a,k y ,n ' 
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V (JV) 

^/„(Qf ))^t ( Xj fcj,)^^, ^ + oQW) + ^I-nWf^L (a, fc,)^ aT (x, fc, + oQW) + h 



We have used 



Y] In{q y )e~ in{u+qyb/2) = e^ [Tj - (tl)+Tj - (u+9! ' b)1 , 



(5.5) 



(5.6) 



n— — oo 



which follows from (2.11). To proceed further, we use the QLA. This amounts to retaining only the gaps that open 
at the Fermi level neglecting those opening above and below the Fermi level. In the QLA, only the term n — N is 
retained in (|5.5|). This leads to 



^t F = -E dxe 



-iaNk-ub 



(5.7) 



In the case of sinusoidal SDWs, Eq. 



]) shows that k y is coupled to k y ± 



ky i 2Qy , 



so that the 



Hamiltonian cannot be straightforwardly diagonalized. For very small 173, the system is not significantly affected by 
umklapp processes (the SDW with wave vector has a very small amplitude), and therefore we expect that the 
spectrum will not be very sensitive to 173. Thus, the most interesting situation arises when umklapp processes are 
strong enough so that the FISDW phase becomes commensurate in the transverse direction (Qy — 7r/6) (which 
implies that N is even). In that case, k y is coupled only to k y + n/b, and both the sinusoidal and helicoidal phases 
are described by the Hamiltonian 



Hmf = (V'+tC^' k v + V 2 ^ $-ii K *i k v ~ */ 2b ) 



K x ,k u 



v F K x 



iNk v b^( N )* 



iNk y b^( N )* 



iNk y b^( N ) _ iNk v bA( N ) 



+t 



+ 1 



y b A 



+t 



+T 
v F K x 



+ 



E 

X 1 ky 



(4> f +l (K x , k y - tt/26), ^(K x , k y + tt/26)) 



v F K x 

_ 6 lNkyb^(N)* 



iNkyb^(N) _ iNkybA( N ) 



e -iNk y b^W* 



-T 



-v F K x 



f $ +r (K x ,k y +Tc/2b) \ 

\^ i {K x ,ky--K/2b) J 



i>+l(Kx,k y - it /2b) 
4>_ r (K x ,ky+v/2b) 



(5.8) 



where K x is now measured with respect to ±k F ± 
I-IbH/v f . In the metallic phase, the dispersion law is 
given by t\ = ±v F K x , e 2 = ±v F K x . The +/— corre- 
sponds to right/left moving electrons and the degener- 
acy is due to spin. In the SDW phase, gaps open at 
the Fermi level, and the dispersion law becomes E\ = 
±[v F Kl + A 1 (fc y ) 2 ] 1 /2 and e 2 = ±[v 2 F K 2 x + A 2 (k y ) 2 ] 1 / 2 
with 



Aife) 2 



A 2 (ky) 2 



\ e -iNk v b^(N) 



lA^ + IA^I 2 



+2\A^y A^y \ cos(2Nk y b + ipi) 



le iNk v b^N) 



. e -WMAW| s 



|AL?I 2 



+2|A (A [ ) A ( ^ ) |cos(2A^fc y 6 + ^2), 



where (pi and <p 2 depend on the phases of 



(pN) 



A. Sinusoidal waves 



(5.9) 



For sinusoidal waves, as follows from (5.9), there are 
interferences between the two SDWs, and the dispersion 
depends on k y . Using |A^' l|/i X ' A ' 
Q4.24|)), we obtain 



\A 2 (k 



V 1 1 nun 



|A^| = |tan(^)| (Eq. 
_ l-\ttm(# N )\ 



|A 1 (fc J/ )| roax |A 2 (fc y )| max l + |tan(^)|' 

For |i?jv| <C 1, which corresponds to a large t/a (section 
U), |A^ I <C |A^p|. The dispersion in the transverse 
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direction is weak. Everywhere on the Fermi surface, the 



gap is of the order of |A^| = |A^|. For larger values 
of |$at|, the dispersion becomes significant. In particular, 
when \$n\ = 7r/4 (which occurs when there is a degen- 
eracy between N and —N, i.e., when £45 = 0), we have 



= IA 



(N), 



and the spectrum becomes gapless.l 



Eqs. (5.9) show that there are 4iV values of k y for which 
Ai or A2 vanish. 

However, the occurrence of the helicoidal phase pre- 
vents the spectrum from becoming gapless. Indeed, the 
stability of the sinusoidal phase requires sin 2 (2$jv) < 
2/3. From (5.10), we then deduce 



|Ai(fc y )| 



V J I mm 



|Al(fc„)|r 



> 0.32. 



(5-11) 



B. Helicoidal waves 



For helicoidal waves, A^j? = A^p = 0. The spec- 
trum is given by E 1 = ±[v 2 F K 2 x + lA^I 2 ] 1 / 2 and E 2 = 



AlT | 2 ] 1//2 . Contrary to the sinusoidal phase, 



±[v%Kl 

the gaps at the Fermi level, A^" 1 and A^\ remain in- 
dependent of the transverse momentum ky. Since the 
sinusoidal phase becomes gapless at \i)n\ = t/4, it is 
natural to expect that, above a certain value of \$n\, the 
system prefers to form helicoidal SDWs in order to lower 
the free energy by opening a large gap on the whole Fermi 
surface. This is precisely the result we have obtained in 
section ^ by considering the Ginzburg-Landau expan- 
sion of the free energy. 



In general, the two gaps |A^| and |A^I ; | differ (they 
are equ al o nly when = or N — 0). As shown in 
section VII, this property gives rise to a kinetic magne- 
toelectric effect. 



AN) l 



VI. QUANTUM HALL EFFECT 

In this section, we use the results of Ref. ^| to study 
the quantum Hall effect. At zero temperature, the offdi- 
agonal conductivity is given by the formula 

CT X y = -ie 2 Y^ I 1TT 7TT {^a\dk v \ll>a) 

-d k Md k M). (6.1) 



The summation is taken over all completely occupied 
bands and the integral is taken over the Brillouin zone. 
{ipaikx, k y )) are the normalized eigenvectors of the Hamil- 
tonian. 

Before calculating a X yin the FISDW phases, we recall 
the main result of Ref. M and generalize it to the case of 
two SDWs. We consider the Hamiltonian 



K*,k„ 



<y 

v F K x Ae"^) \ ( j, + (K x ,ky) , .. 
Ae^W -v f kJ \ ^{K x ,k y ) ' [b - Z) 



describing spinless electrons in presence of a charge- 
density wave at wave vector (2kp, Q y = 0). Here A and 
4>(k y ) are real. K x = k x =p k F for ± electrons. The 
spectrum is given by E+ — (v F K x + A 2 ) 1 / 2 for the up- 
per (empty) band and £L = —{v F K x + A 2 ) 1 / 2 for the 
lower (filled) band. The wave functions are defined on 
the Brillouin zone torus \k x \ < k F and < k y < lixjb. 

Consider now a wave function \ip(k x °\ k y )) in the lower 
band. When k x is changed along the closed line encir- 
cling the torus at fixed k y , the wave function changes and 

transforms intp,,e^' Cy ' ) \ip(k x °\ k y )) when we return to the 



sta rtin g point .cn£3 Noting that the second term in the rhs 
of ( |6.l| ) vanishes (since the Hamiltonian is the same for 
ky = and k y — 2ir), we obtains 



>xy 



An 2 



[0(2tt/6)-^(O)]. 



(6.3) 



We now consider electrons with spin in presence of two 
SDWs (with zero transverse wave vectors: Q y = 0): 



Hmf= ^2 vP+\( K x> k y)> $-i( K x,k y 



v F K x A ie -*( fc *) \ / i> +1 (K x , k y ) 
A ie *( fe «) -v F K x ) \ i,_ L (K x ,k y ) 

( v F K x A 2 e"^( fc «) \ / i> +i (K x , k y ) 
y A 2e *( fe «) -v F K x ) \ ^(K x ,k y ) 

(6.4) 

where Ai, A 2 , <j>i, and 2 are real. We apply the same 
reasoning as before to the wave function \ip^(k x °\ k y )) 
of the lower band. After the trip around the Brillouin 
zone, we obtain the wave function e 1 ^ 1 ^^ i{k x °\ k y )). 
Repeating this procedure once more, we obtain the wave 
function e <<M*v | ^ T (fc|?\ fc v )). This yields 



<r*v = 7- [</>! (2tt/6) + fa (2tt/&) - 0i (0) - 2 (0)] . (6.5) 



It was shown in Ref. |9| that ( |6.5| ) does not depend on the 
value of the transverse wave vectors of the SDWs. More- 
over, the global contribution oLthe gaps (if any) located 
below the Fermi level vanishes!] We are therefore now in 
a position to calculate a xy for the FISDW phases, using 
the Hamiltonian derived in the QLA (section [v]). 
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A. Sinusoidal waves 



Comparing ( |5.8| ) and (|6.4|), we deduce 

= -e^^lA^I - e- mk * b \AW\. (6.6) 



We have used |A 



OJV)| 
+ T 



|A_| |. When fc y varies from 



to 2tt/&, the change ^i(2tt/&) - ^i(O) = 2 (2tt/&) - </> 2 (0) 
is determined by the term with the largest amplitude in 

©. Since lA^/A+i?! = |t| < 0.518 in the sinusoidal 
phase, we obtain 



0i(2tt/6) + 2 (2tt/6) - 0i (0) - 2 (O) = -4ttA. (6.7) 
This yields 



e 2 



(6.8) 



where we have restored the dimensional Planck constant 
h. Thus, in the sinusoidal phase, the QHE is determined 
by the SDW with the largest amplitude. 



B. Helicoidal waves 



For helicoidal waves, we deduce from (5.7) 



/^ ie -i<t>i(K) = _ e WMAW, 



(6.9) 



This yields 0i(2tt/&) ~ <M0) = ~[M^/ b ) ~ <M°)] so 
that 



0. 



(6.10) 



Thus, we come to the conclusion that the QHE vanishes 
in the helicoidal phase. 



VII. MAGNETOELECTRIC EFFECT 

A magnetoelectric effect may exist if time-reversal and 
space-inversion symmetries are broken.Ea Gor'kov and 
Sokol found the kinetic magnetoelectric effect for a single 
helicoidal SDW. £3 The effect also exists in the presence 
of two helicoidal SDWs of opposite chiralities, provided 
their amplitudes are not equal. An electric current j x 
along the chains induces a uniform magnetization <5M 
along the vector n that characterizes the spin polariza- 
tion of the helicoidal SDWs. In our case, the vector n is 
parallel to the magnetic field H, which is oriented along 
the z axis, thus 



SM Z 



(7.1) 



Here 5M Z is the additional spin-magnetization density 
induced by j x in excess of the magnetization density M z 
induced by the magnetic field without j x . The effect can 
be understood by considering the spectrum of electronic 
excitations in the helicoidal FISDW phase shown in Fig. 
^. The +kp electrons with spin up and the — kp electrons 
with spin down have the energy gap |Ajv.+ |, whereas the 
electrons with spin down and the — kp electrons with 
spin up have the different energy gap | A^ _ |. To produce 
a current j x along the chains, electrons need to be trans- 
ferred from — kp to +kp (we assume that the electric 
field is weak enough so that the SDWs remain pinned by 
impurities). For |Aat :+ | ^ |A^_| (|Ajv,+ | ^ |Ajv,-I if 
<4t 7^ and N ^ 0), this redistribution of electrons af- 
fects up and down spins in a different way, which results 
in a uniform magnetization M z . Denoting the deviation 
of the distribution function of electrons with spin a and 
momenta near akp from the equilibrium one by Sn ajlT , 
we have 

j x = ev F ^(<Sn +jCr - Sn- ia ), 
5M : 



5Ms 
2 



^(5n a ,t -Sn a ,i). 

a=± 



(7.2) 



Here we denote the electron gyromagnetic factor by g 
(assumed to be equal to 2 in the preceding sections). At 
low temperature, T <C |Ajv,+ |, |A_y |, the electrons are 
excited solely above the lowest energy gap (Ajy _ in Fig. 
^). This implies that Sn + ^ ~ — and 



SM Z 

jx 



9VB 

2evF 



Eq. (7.3) can be rewritten as 

SM Z gl 



Ms 



2evpL v L z 



(7.3) 



(7.4) 



where / = j x L y L z is the current passing through the 
sample of cross-section L y L z . For / ~ 1/iA (which is 
slightly below the critical current for the depinning of 



the SDWsEI), L y L z ~ 1 mm 



2, we obtain 
\6M Z 



(IB 



2 x 10" 



vp ~ 3 x 10 m/s, and 



A' 3 . (7.5) 



This should be compared with the ground-state magne- 
tization density 



\M Z \ 



10~ 8 A' 



(7.6) 



that we obtain from Ref. |7[ Thus, we obtain 
\5M Z \/\M Z \ - 2 x 10~ 7 . 

The reentrant phases N — are somehow special 
since |A^j| = |A^°j| independent of t±b- Consequently, 
these phases do not exhibit the magnetoelectric effect, 
although they are helicoidal. 
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VIII. COEXISTENCE BETWEEN SUCCESSIVE 
PHASES 

It has been shown by Lebed' that under certain circum- 
stances, umklapp processes can lead to the simultaneous 
existence of two successive sinusoidal phases N + (3 and 
N .113 The system then evolves from the phase N + (3 to 
the phase N via a region of the phase diagram where 
both the phases N + (3 and N exist. The transitions to 



J 



the coexistence region are of second order. 

In this section, we reconsider the problem of the coex- 
istence between two successive SDW phases, distinguish- 
ing between sinusoidal and helicoidal waves. We analyze 
all possible cases (depending on the polarization of the 
SDWs). 

The coexistence of two phases implies the si mul tane- 
ous fo rmation of four SDWs. Therefore Eqs. (4.7) and 
(4.10) should be replaced by 



A U„ „, ^ — V V X (n A(P( N +"f)) -L n A ( - P( - N+ ^h P - la Q i v (N + l))b (h+h)/2 T / Jn( N +"f)\ 

/\ aa (Wi,W 2 ) — } J 2 , °k 2!c .k l!c +aQ { I {N + ' ,)) ^ g ' 2 ^ +93^-aa y I a(h-l 2 ){Qy )■, 

p=±7=0./3 

V V 8 , , w+ »8, , fw , , e -«QS, P(JV+ "' )) b[2il-ap(JV+ 7 )]/2X(p(JV+ 7 )) 

Z / / / fe 2a: ,fe la: +aQi P<N+ "' ))0 '2,h-ap(W+7) e ^aa 
p=±7=0,/3 

The free energy becomes 

pint 



A atr (wi, w 2 ) 



QLA 



(8.1) 



F 



N,N+0 



= F, 



N 



F, 



N+p 



F 



N,N+0i 



(8.2) 



with the interacting part 

2KJ2{\^\ 2 



pint 



|A 



a] I 



IA 



a] I 



IA 



(JV+/3)|2 
qT I 



-2^cos[(iV + 0)Qi N) b - NQl»+»b] ^(aWAS'A^A™* + c.c.) 



^E[ 2 E l A if MM* +cos(/3Q(»)6)[(Ai°)) 2 + (Ag) 2 ]Ag)*Ai' 



(0h 



+ c.c. 



(8.3) 



Without loss of generality, we have assumed that the phase N + /3 is not the last phase of the cascade (i.e., we do not 
consider the case N + (3 = and Q y N+ ^ — n/b). In order to study the possible coexistence of phases N and N + [3, 
we consider a region in the phase diagram where the phase N is more stable than the phase N + (3 (Tc N ^ > Tc N+ ^). 
Assuming that the order parameters of the phase N + [3 are infinitesimal, we derive an effective free energy for this 
phase, from which we conclude about the coexistence of the two phases. The same reasoning is applied to regions of 
the phase diagram where Tc < Tc N+ . 



A. Sinusoidal waves 



1. N / 



Introducing the order parameters u« , Va , u\^ + ^\ and Va^^ p> defined in section ttVAl, and setting Va" ' — 



(N+f3) 



,w _ 



0, we rewrite (EM (for TV ^ 0) as 



pint 

r N,N+f3 



2K^[cos 2 (^)|uW| 2 +sin 2 (^)|4 A ' ) | 2 ] [cos 2 (^ +/3 )| M ^+«| 2 + sin 2 (^ +/3 )|uf + Y 

a 

+Kcos[(N + 0)QW b - NQ^+f%] sin(2^) sm(2d N+0 )(uf ] u W u^ +fl) *u {N+fl) * + c.c). (8.4) 



\u W \ and \ui N+f3) \ 



\u_ + \. We consider a region of the phase diagram 



For sinusoidal waves, we have \u 

where T<P^ > Tc N+ ^ and assume that u^ + ^ is infinitesimal. To lowest order in til"^', the effective free energy 
for the phase N + (3 is 



TTlCff 



= 2\\ N+P) \u { ? +P) \ 2 +4#|uf >|>f +P) \ 2 



+Kcos[(N + P)Q ( y N) b ~ NQ y N+ ^b] S m(2ti N ) sm(2V N+p )(u^ } u i _ N) u ( f +m *u (N+m * + c.c.) 



(8.5) 
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where is not changed by the infinitesimal u K a +p> and is therefore given by (pL2q) for T < T^ ' ■ The free energy 



XN+0) 



is minimal if the phase of ui N+P) is such that 

i(N + p)Q<f) h - NQ y N+ ^b] sm(2$ N ) sm(2$ N+ p)u^ui N \i N+f3 K { _ N+f3) * 
[(N + 0)QWb - NQW+Vb] sin(2^) sin(2^ +/3 )| \u^u^ + «| 2 , 



cos 



(8.6) 



which yields 



Ff +0 = 2\u + 



(JV+/3).2 



X (N+P) + 2K \u { f } \ 2 - K\uf } | 2 | cos[(iV + 0)QWb- NQ y N+ ^b] sm(2V N ) sin(2t? A r+ / a) |1 . (8.7) 



A second order phase transition to a phase where both 
and u+ + are nonzero occurs if the coefficient of 
| u (7V+/3)| 2 - n becomes negative. Since in all cases, 

Qy , Qy + ~ 7r/& (see section R), we make the approx- 
imation | cos[(N + (3)Q ( v N) b-NQY +l3) b}\ ~ 1. Using (see 
appendix |(]) 



we can write the coexistence condition for T < Tc N ^ as 



T-T, 



N+/3) 



T, 



(N+0) 



< 



T-T, 



(N) 



T 



(N) 



-x+, 



(8.9) 



where 



|sin(2^ JV )sin(2^ JV+(3 )| 



1 



: sin 2 (2tf 



N) 



where x = X-, x+. Since T — Tc N ^ < 0, the coexistence is 
possible only if X- < 1 or x + < 1. The existence of lin- 
early polarized SDWs requires sin 2 (2i?Ar),sin 2 (2^Ar + ^) < 
2/3 and therefore implies a?— Imin = z+lmin = 1< Con- 
sequently, there is no coexistence between phases for 
sinusoidal j-ajaves. Thus, our results invalidate Lebed's 
conclusion!!! concerning the coexistence of two successive 
sinusoidal phases. 

If we discard the existence of helicoidal waves, then 
the coexistence between phases may be possible depend- 
ing on the geometry of the Fermi surface. Let us take 

for instance Q y = Qy + — n/b and t^b = 0, which is 
the case considered by Lebed'. Then -d^ = $N+p = 7r/4. 
The region of coexistence does exist and is determined 
by 



3T — Tc 



(N+/3) 



2 j,(N+0) 



< 



T-T 



(N) 



T 



(A r ) 



2T -T x c 
< -- 



T 



JV+/3) 



.14) 



(8.10) This is precisely the result obtained by Lebed' 



.0 



For the coexistence in regions of the phase diagram 
where Tc N ^ < Tc + , we obtain by the same method 
the condition (for T < T^ N+P) ) 



T-T, 



(N) 



< 



T-T 



(N+f3) 



T 



(N+0) 



where 



|sin(2T?jy)sin(2i?jv +/3 )| 
l + isin 2 (2^ +/3 ) 



(8.11) 



.12) 



$n, $N+p are determined by An, and B, and 
therefore depend on temperature and magnetic field. 
To verify tha t the regio n of the phase diagram defined 



by Eqs. ( |S.9|) and ( 3.11 ) exists, we consider a point in 



the (H, T)-plane slightly below the intersection point 



of the two curves T, 



(N) 



and T, 



(N+/3) 



T 



(A') 



For this point 



rp(N+p) an( j ^ e coex j s ^ ence conditions can be 



rewritten as 



x) 



T-T, 



(JV) 



T 



(N) 



< o, 



(8.13) 



2. N = 0, 



l(0) 



= Tl/b 



We can study the coexistence of the last two phases 
of the cascade f ollowi ng the same procedure. Using the 
result of section IV B, we have 



+2/sTcos(/37r) sm{2d i3){u 2 uf ] *u m * + ex.). (8.15) 

Choosing the phases of and u_ in order to minimize 
the free energy, we obtain 

F^=4if[2-|sin(2^)|]| M | 2 | U f | 2 . (8.16) 
We deduce that the coexistence region is determined by 



1 T-T 



(0) 



x^ 
where 



< 



T-T 



(o) 



T 



(0) 



<x'_ 



T-T 

rr(0) 



(0) 



.17) 



c , _ 2-|sin(2^)| 

l + isin 2 (2^)' 
c' + =2 - |sin(2^)|. 



(8.18) 
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The coexistence occurs if x', < 1 or x'_ < 1. Since 
sin (2i?0) < 2/3, x'_| m i n ~ 0.89. In principle, the coexis- 
tence between the phases N = and N — 0, Q y — n/b 
is therefore possible in the region of the phase diagram 



where T^ v > T^> provided that x'_ < 1. If (3 = 1, then 
our numerical calculations show that sin (2$^) <C 2/3 
(the latter inequality holds for all phases with N odd). 
Therefore, the coexistence region does not exist if the 
phase N = is preceded by the phase (3 — 1. In the case 
where the phase N = is preceded by the phase (3 — 2, 
the coexistence region exists provided that x'_ < 1. This 
situation however requires strong umklapp scattering and 
is th erefo re quite unlikely. (Note that Jar sin 2 {2'dp) = 1, 



Eq. (8.17) agrees with Lebed's results. 113) 



AN+0) 



0, this would not change the result). Eq. (|8.3|) 



then yields 



pint _ K '| W W„(^+«|2 



.22) 



In the region Tc > Tc N+ , we find the effective free 
energy for the phase N + (3 (for T < T c (Ar) ) 



pelt _ I 
^JV+/3 — l M + I 



J. _ T (N+P) rp _ T (N) 



T (N+P) " T (N) 



(8.23) 



5. JV + (3 = 0, ?4 n/b 



The coexistence between the phase TV and the reen 
trant phase N + (3 = 0, Q y N+ ^ ^ n/b can be stud- 
ied by setting (3 = —N. This case is special since the 
phase N + (3 = 



0, Qy N+ ^ j^z n/b is always helicoidal 



when g 3 ^ (see section |v|). Using 



lu^l and 



,(°) 



=0, we have 

rmt _ O!/l.,(0)„mi2 



FT = 2K\u^u\ 



.19) 



In the region T { C N) < T c w , wc find the effective free en- 
ergy for the phase N (for T < T c (0) ) 



,(0) 



The coefficient of 



,(iV+/3)|2 



is always positive so that 
there is no phase coexistence. In the region T c < 
we find the effective free energy (for T < 



(N+/3) 
(N+0) 



T-T, 



(AO 



T-T, 



{N+fi) 



T 



(AT) 



JV+/3) 



which shows that there is no coexistence. 



(8.24) 



i^f = | u f>| 2 jV(0) 



T —T, 



(N) 



T — T, 



(0) 



r, 



(0) 



.20) 



The coefficient of lu^l 2 is always positive so that there 
is no phase coexistence. In the region < Tc N \ we 



find the effective free energy (for T < T< 



f; 



■off 



,(o) |2 JV(0) 



T-T, 



(0) 



(o) 



T-T, 



(N) 



l + ±sin 2 (2iV) T W 



which shows that there is no coexistence. 



.21) 



g. TV = 0, Q y = tt/o 



For the coexistence between the helicoidal phase N = 
(3 and the sinusoidal phase N — 0, Q y — n/b terminating 
the cascade, we have 



0,0 



This yields (for T < T c (0) ) 



Ff 



,03)12^(0) 



r-T, 



(/3) 



T-T, 



(0) 



(0) 



(8.25) 



(8.26) 



B. Helicoidal waves 

1. N / 0, or N = and qJ, 0) / 7r/b 

First we study the coexistence between the two heli- 
coidal phases N and N + (3. For simplicity, we consider 
only the case i?jv = $n+/3 = tt/4. For helicoidal waves, 
we have u_ = and u^ + ^' = (we could also choose 



in the region T c (0) > T c (/3) , and (for T < T c w ) 



(0h 



T-T, 



(0) 



T-T, 



(0) 



T 



(0) 



(/3) 



.27) 



in the region T^ " 1 < Tc^ . The coexistence is therefore 
not possible. 
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IX. SUMMARY AND CONCLUSION 

Our main results can be summarized as follows: 
In presence of umklapp processes, the instability of 
the metallic phase at the temperature Tc N ^ corresponds 
to the formation of two SDWs, with quantized longi- 
tudinal wave vectors Q x = 2k F + NG and Q x = 
2kp — NG. For very weak umklapp scattering, both 
SDWs are incommensurate in the transverse direction 
{Qy N ^ — —Qy N ^ 7^ except when N = 0. If we label 
each phase by the integer TV corresponding to the SDW 
with the largest amplitude, we have sgn(./V)=sgn(i2b). 

The amplitude of the SDW at wave vector Q x is van- 
ishingly small. The quantum Hall conductivity is de- 
termined by the SDW with the largest amplitude, i.e., 
a xy = -2Ne 2 /h. 

For even N, there exists a critical value of <?3 (typi- 
cally (73/32 of the order of a few percent) above which 
the system prefers to form two transversely commensu- 
rate SDWs (Q y N) = -Q y N) = n/b). When t ib > 0, 
the SDW with the largest amplitude is then determined 
by sgn(iV)=— sgn(t2h): the QHE changes sign. The 
two SDWs have comparable amplitudes when is not 
too large (£4;, < 1.5 K). When the umklapp scattering 
strength increases, the first negative phase to appear is 
the phase N = — 2. 

Umklapp scattering also tends to suppress the phases 
with an odd N and produce some reentrances of the 
phase N — within the cascade. Unlike the last phase of 
the cascade, the reentrant phases N — are incommen- 
surate (Q y 7^ n/b). 

The negative phases are likely to become helicoidal 
when the umklapp scattering strength is further in- 
creased. Experimentally, this situation could be achieved 



by decreasing pressure. The appearance of these heli- 
coidal phases is entirely controlled by . The QHE van- 
ishes in the helicoidal phases, but a magnetoelectric effect 
appears. These two characteristic properties may be uti- 
lized to detect the helicoidal phases experimentally. The 
reentrant phases N — are always helicoidal but do not 
exhibit the magnetoelectric effect. 

In the sinusoidal phases, umklapp processes modulate 
the gap on the Fermi surface as a function of k y . When 
|A(/cy)| m i, 1 /|A(fc !y )| max ~ 0.32, the sinusoidal phase be- 
comes unstable against the formation of a helicoidal 
phase. ■— ■ 

The conclusion of Lebed'E-D that, in the presence of 
umklapp scattering, adjacent FISDW phases are sepa- 
rated by two second-order phase transitions and an in- 
termediate phase with coexistence of four SDWs is incor- 
rect, because he did not consider helicoidal SDWs. 

In conclusion, the consideration of umklapp scat- 
tering naturally explains the appearance of negative 
FISDW phases in quasi-one-dimensional organic conduc- 
tors. These phases are characterized not only by a sign 
reversal of the QHE, but also but the simultaneous pres- 
ence of two SDWs with comparable amplitudes (provided 
that is not too large). This leads to the possible sta- 
bilization of helicoidal phases. But even in the sinusoidal 
phases (which are the ones that have been observed up 
to now), we expect the presence of two SDWs to give rise 
to new physical properties. 
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APPENDIX A: 



We calculate in this appendix the pairing amplitudes 



A a<7 (w ly w 2 ) = J dVC;(r)C 2 (r)A a ,(r) 

P=± J 

Using -4k F + Q x pN) = -Q x pN) and Q y pN) = -Q ( y pN \ we have 

P =± J 



Q P N r 



Using (2.11), we have 

x •* m 

=^, fclI . +QQ ->^ MQHb(/i+i2)/2 ^(^M(^r ) ) 



a -iaQ i ;f N ' } bm 



H2—m 



e 



Since I n (q y ) = I n (-q y ), (|A3|) leads to (J4J 



(Al) 



(A2) 



(A3) 



17 



APPENDIX B: 



In this a ppen dix, we give the main steps in the calculation of the fourth order term of the free energy leading to 
( 4.13 ) and ( |8.3| ). We consider the case where there are simultaneously four SDWs at wave vectors Qn, Qn, Qjv+/3, 

Qn+/3 : 



F 



(4) 



K 



N 



2 E E E 



x exp 



X (pi(JV+7i)) X (P2(N+ l2 ))* X (p 3 (N+ l3 )) X (p 4 (AT+ 7 4))* 
at at a| at 



{»(6/2) [pi(JV + 7i)Q^ l(JV+7l)) ~ [2pi(JV + 71) - P2(N + l2 )]Q^ 2(N+l2)) 
+ [2 Pl (N + 7l ) - 2 P2 (N + 72 ) + p 3 (iV + 73 )]Q^ 3(Ar+73)) - P 4 (^V + 74 )Q< / P4(JV+74)) ] } 



x5 



(Pl-P2+P3-P4)JV+pi7i-p 2 72+P373-P474 



>° E *Q<« 



(JV + Tl)) _q(P2(- n + T2)) _|_q(P3( w +T3)) _Q< j P4 ( N + T4)) iT j27r/6 ' 

(Bl) 



We write F 



(4) _ F (4) 



AT 



JV 14 



(4) 



JV 13 



Fn^+Fn^i- F n'U corresponds to the case where all the (pj, 7 j) in (|Bl|) are 



identical (p x = p 2 = • • •, 7l = 72 = • • •) 



■fjv la corresponds to the case where three of the (pj, 7 j) are equal and 



different from the fourth one... F^ |i corresponds to the case where all the (pi, 7 i) are different. Skipping the details 



of the calculation, we only give the final result for F 



(4) I 

N !*• 



J- ) 



AT 13 

-(4) I 
JV |2 



AT |1 



= - y ia^' 



(p(AT+ T ))|4 



a,p,7 



<w(E ( W„ 



ky [a^ia^pas 



^E E iAir +7))A S 

a py^p'-y' 



(o)|^(B).2X(0> 

ay I at 1 

a (p(AT+ 7 )) r (p'(Ar+ 7 '))i2 



A (5 M ( °T ) | 2 A ( t 



cos(/3Q(°>6)[(A^ + (A^]A^*A^* + 



>(°),2 



2 y-| A ( P /3)|2 A (0) A (0)* 
/ ; I a| I ay at 
P 

(0)n2i a(/3)*a(W* 



1 ,7 

2 



c.c. 



(A<°>A<°>) 2 

v a| a>\ / 



C.C. 



2^cos[(7V + /3)Q™& - tfQf +«6] ^(Aj'Aj'Af^A'^^ 



a? 



(B2) 



APPENDIX C: 

In this appendix, we calculate \[ N ^ for T ~ Tc N \ Since A x vanishes for T = Tc N \ we have 



Af'~(T-Tf)) 



9Ai 



(AT) 



<9T 



(T-T C W) 



(Cl) 



From (2.19) we deduce dB/dT = and 



An — An\ t (n) 



T-T C W d X f(Q N ) 



T 2 



dT 



(C2) 
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with 



dx ( +\^N) = Y / l^Q^)^- {2kF + {N ~ n)G) 



dT 



dT 



(C3) 



Eq. ( gJO[) yields 



d X ^(2k F -nG) _ N(0) j 



dT 



Re 



inujr , ( 1 inuj r \ 



4ttT 2 V2 4ttT 



2 I T 

where is the derivative of the digamma function. Using $>'(z) ~ 1/z + l/2z 2 for |z| 3> 1, we obtain for lu c ^S> T 

d X 1 ^(2k F ~nG) N(0) 



(C4) 



<9T 



2T 



■4,o + 0(TK). 



To leading order in T/u> c , we therefore have 



3x1 0) (QaO _ r2 JV(0) 



<9T 



N 2T ' 



from which we deduce 



In the same way, we can show that 



This yields 



A N = A N \ Tr +^(T~T^) 



2T, 



(C5) 



(C6) 



(C7) 



(C8) 



(C9) 
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FIG. [j]. Diagrammatic representation of the integral equation for the susceptibilities x++ an d Xh m the RPA. The solid 

(dashed) lines represent electrons on the right (left) sheet of the Fermi surface. The arrows indicate the spin direction. The 
wavy lines correspond to forward (52) or umklapp scattering (ga). 




H(T) 

FIG. ^| Phase diagram for r 



(<72 — 0.38 and 173 = 0) and t^, = 0.75 K. Top picture: Transition temperature T c 



(JV) 



(solid line). The dashed line gives the value of the integer N. The QHE is determined by a xy — —2Ne /h. Bottom picture: 
Transverse wave vector Q y maximizing the transition temperature ■ 
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H (T) 

FIG. |. Phase diagram for r = 0.025 (g 2 ~ 0.37 and p 3 ~ 0.01) and t 4b = 0.75 K. The bottom picture shows sm 2 (2-& N ) 
(solid line), \"f\ (dashed line), and I7I (dotted line). sin 2 (2 - i?jv) determines the polarization of the SDWs and I7U7I give the 
ratio of the amplitudes of the two SDWs (see text for a precise definition). (These quantities are not shown in the last phase 
of the cascade (iV = 0, Q y — ir/b) which contains a single circularly polarized SDW.) A negative phase (iV = —2) appears 
in the cascade. The two SDWs have comparable amplitudes in that phase: I7I, I7I ~ 0.3. All the phases are sinusoidal (since 
sin 2 (2t?jv) < 2/3) and the Hall effect is quantized: a xy = —2Ne 2 /h. 
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FIG. ^. Phase diagram for r = 0.025 (c/2 — 0.37 and 33 ~ 0.01) and tn> = —0.75 K. All the phases are positive, but the phase 
N — 2 has split into two subphases, one of which is commensurate in the transverse direction (Q y = ir/b). 
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FIG. ^ Phase diagram for r = 0.03 and = 0.3 K (the horizontal dot-dashed line corresponds to 2/3). When t^, is reduced 
from 0.75 to 0.3 K, the negative phase N = —2 becomes helicoidal (sin 2 (2i9_2) > 2/3). The QHE vanishes in that phase, but 
a magnetoelectric effect appears. 
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FIG. []. Phase diagram for r = 0.06 and = 0.75 K. Two negative phases, N = —2 and N = —4, can be observed. The 
phase iV = —2 has split into two subphases: one is helicoidal (sin 2 (2tf_ 2 ) > 2/3), one is sinusoidal (sin 2 (2tf_ 2 ) < 2/3). 
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FIG. [?]. Phase diagram for r = 0.2 (§2 — 0.32 and 33 ~ 0.06) and = 0.75 K. Only phases with negative even N survive 
when umklapp scattering is strong enough. Some reentrances of the phase N = (with Q y 7^ n/b) also appear within the 
cascade. 
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FIG. g. Free energy Fn as a function of |u+ | and \u_ | for different values of the parameter x = sin 2 (2i?jv). When x < 2/3, 
the minimum of -FV corresponds to (sinusoidal phase). When x > 2/3, there are two minima located on the 

line = and u_ = (helicoidal phase). When x = 2/3, the minima are infinitively degenerate. 
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FIG. y. Spectrum of electronic excitations in the helicoidal FISDW phase. The solid (dashed) line corresponds to up (down) 
spins. Only the gaps at the Fermi level are shown. For clarity, we have not shown the Zeeman splitting. 
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